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Abstract 

Evidence for two types of relatively large amplitude MKD waves upstream 

and downstream of quasi-parallel forward and reverse interplanetary shocks t* 

presented. The first mode is an Alfven wave with frequencies (in the 

spacecraft frame) in the range of 0.025 - 0.07 Hzi. This is a left-hand 

polarized mode and propag'Wes within a few degrees of the ambient magnetic 

field. The second is a fast MHD mode with frequencies in the range of 0.025 - 

0.17 Hz, right-hand polarization and propagating along the magnetic field. 

These waves are detected principally in association with quasi-parallel shocks 

( 9 < 45°). The Alfven waves are found to have plasma rest frame frequencies 

on 

in the range of 1.1 - 5.3 mHz with wavelengths in the order of 4.8 x 10 8 - 2.7 

x 10* cm. Similarly, the fast MHD modes have rest frame frequencies in the 

range 1.6 - 26 mHz with typical wavelengths about 2.19 x 10* - 9.51 x 10* cm> 

The magnetic field power spectrum in the vicinity of these interplanetav 
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shocks is much steeper than f at high frequencies. The observed spectra 

-2.5 -4 

have a high frequency dependence of f to f . A peculiar feature of the 
fast mode identification in one event is the large correlation observed be- 
tween |B| and proton density p for field aligned propagation. This appears to 
be a nonlinear effect, second order in the wave amplitude. An interpretation 
of these observations is given in terms of the electromagnetic ion beam insta- 
bility. Both resonant and nonresonant interactions need to be considered to 
account for the polarization and spectral content of the observed fluctuations. 
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1. Introduction 

Recent experimental observations and theoretical modeling have pointed out 
the Importance of plasma waves upstream of interplanetary shocks in under- 
standing the structure of collisionless shocks and the origin of energetic 
particles often observed in association with shocks [Tsurutani et al., 1983; 
Acuna et al., 1981). These observations have detected two types of upstream 
waves. A high frequency whistler mode with frequencies ag measured in the 
spacecraft frame between 0.2 and 2 Hz and a low frequency fast MHD mode near 
50 mHz. Both wave modes are observed to have circular or elliptical right- 
hand polarization (in the plasma frame), and propagate within 15° of the mean 
magnetic field direction. Similar observations of wave phenomena upstream Of 
planetary bow shocks, including observations of left-hand waves, have been 
discussed by Barnes [1970], Fairfield [1969, 1974], Russell et al. [1971], 
Hoppe and Russell [1983], Smith et al. [1983], Goldstein et al. [1983] and 
Smith et al . [1984] . 

In this paper, we present preliminary results of an investigation of 
magnetic fluctuations seen upstream of two interplanetary shocks. The 
spectral analysis includes calculation of the normalized reduced magnetic 
helicity spectrum cr^Ck), the normalized reduced cross helicity spectrum a c (k), 
and the Alfven ratio r (k) as discussed by Matthaeus and Goldstein [1982]. 

A 

Minimum variance methods are used to compute wave polarization as a function 
of frequency. The Taylor "frozen- in- flow" hypothesis is assumed [Taylor, 
1935, 1938; Matthaeus and Goldstein, 1982] to convert frequencies to wave- 
vectors. Some of the basic properties of the waves including the probable 
mode of propagation in association with both quasi-parallel forward and 
reverse shocks are described. A comparison with previous results on the 
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generation of waves at interplanetary and planetary shocks is presented. 

Section 2 contains a discussion of the spectral techniques. Section 3 
contains a discussion of the linear theory of the electromagnetic ion beam 
instability which we utilize to explain the excitation of the observed 
fluctuations. The results are discussed and summarized in section 4. 

2. Observations 

The Voyager (Gf> FC) magnetometer and (MIT) plasma teams have compiled a 
list of interplanetary shocks from launch (October 1977) to 1980 for time 
intervals when simultaneous observations are available. The plasma parameters 
were determined from either moment calculations or from a Gaussian-fit 
procedure as discussed by Bridge et al. [19771. We used only the moment 
calculations which were computed from data sampled every 12 s during the time 
periods discussed in this paper. The magnetometer aboard Voyagers 1 and 2 has 
a much higher time resolution [see Behannon et al., 1977]. Magnetic power 
spectra and magnetic heliclty spectra were computed using 1.92 second average 
magnetic field data. Consequently, these spectra extend to higher frequencies 
than the plasma data. Shock normals were calculated from the plasma-magnetic 
field data using a single spacecraft method of shock normal estimation develop- 
ed by Lepping and Argentiero [1971] and improved by Acuna. The time periods 
of interest in this paper are 0800 - 1000 UT on January 29 and 1850 - 2000 UT 
on February 3, 1978. 

January 29, 1978 — Upstream 

Figure 1 shows a plot of the magnitude and components of the magnetic 
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field and plasm* bulk velocity for January 29 in Che RTN (radial, tangential 
and normal) heliocentric coordinate system together with the proton number 
density and temperature. Note the presence of a quasi-parallel reverse shock 
at about 0918:17 UT. The shock normal components for this event are n*(-0.92, 
-0.35, -0.17) and the angle between the normal and the average upstream 
magnetic field is about 11°. The shock structure and its properties for this 
event has been recently investigated by Scudder et al. [19831. For the 
spectral analysis calculation we have selected the upstream and downstream 
regions (excluding the shock itself) which correspond to the subintervals from 
0919:10 to 1000 UT and 0800:45 to 0914:39 UT, respectively. 

For the upstream subinterval the components of the average magnetic field 
are (-1.13, -0.56, -0.46) nT, and the average sr,ilar wind speed is 366 km/s. 
The fluctuations in all the components and in the magnitude are rather large, 
typically ( SB) ms / K B 0 ' > I “ 0*40. As usual when data is available from only a 
single spacecraft, only reduced (one dimensional) spectra can be determined 
[Batchelor, 1970]. We use the fast Fourier transform technique with 26 
degrees of freedom to compute these spectra. Details of our analysis techni- 
ques can be found in Matthaeus and Goldstein [1982]. In Figure 2 we show the 
magnetic field power spectrum corresponding to this subinterval. Note that it 

posseses an f ^ power- law dependence at high frequencies, which differs 
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substantially from the typical ambient solar wind behavior of f (Joklpii 

and Coleman, 1968; Matthaeus and Goldstein, 1982]. The bulge in the spectrum 

centered about f - 8 x 10 ‘ Hz reflects the presence of quaai-monochromatic 

fluctuations in the upstream region of this shock. 

The normalized reduced magnetic helicity spectrum a ffl (k) = \ kj H m (k) /E b (k) 

is plotted in Figure 3a. F.,(k) and H (k) are the reduced spectra of magnetic 

d m 

energy and magnetic helicity, respectively. E^(k) is the trace of the 
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v* 

spectral matrix (k) defined from [cf. Matthaeus and Goldstein, 1982]: 

s i j 1 *(k 1 ) - / dkjdkj S ij (k l ,k*,k, ) 

The magnetic helicity is defined by 

H = /d 3 x X*S - /d J k H (k) 
m m 

The reduced helicity spectrum H (k) is computed from [Matthaeus and Goldstein, 
1982]: 

H m (k) - 2 Im S„ r (k)/k 

where the components "23" correspond to "TN" . The total magnetic helicity H 

ra 

* /dk H (k), is a measure of the lack of mirror symmetry of the magnetic 
field and determines its topological handedness or "knottedness" [Moffatt, 
1978; Matthaeus and Goldstein, 1982]. The magnetic helicity spectrum will be 
positive for left-hand topological structures and negative for right-hand 
topological structures. The sign of the magnetic helicity is directly related 
to the sense of polarization as defined in optics [Jackson, 1962]. A discus- 
sion of the relationship between the magnetic helicity and the polarization 
can be found in Smith et al. [1983], An important fundamental property of the 
magnetic helicity is th&t because it is a Galilean invariant, once determined 
in the spacecraft frame it is also known in the plasma frame (subject to the 
constraint that one is dealing with MHD phenomena) . Therefore the sense of 
polarization determined from it corresponds to the polarization in the plasma 
frame. Note in Figure 3a that the magnetic helicity increases with frequency 
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becoming positive (o m ~ 0.8) for frequencies in the range of 0.025 - 0.07 Hz, 
indicating the presence of left-hand helices in the magnetic field. This 
magnetic helicity spectrum is distinctly different from the randomly ocillat- 
ing spectra normally observed in the solar wind far from planetary and inter- 
planetary shocks [Matthaeus and Goldstein, 1982; Smith et al., 1983). 

An eigenvalue (minimun variance) analysis of this interval was performed 
in which the spectral matrix was rotated into an eigenvalue Coordinate system 
at each mode. This yields the degree of polarization and ellipticity of the 
fluctuations as a function of frequency. The smallest eigenvalue is associat- 
ed with the direction of minimum variance of the fluctuations for each Fourier 
mode. This is the direction of ±1c under the, assumption that decomposition 
into plane waves is appropriate. 

Figure 3b-d shows the results of this analysis. In Figure 3b we plot the' 
degree of polarization D, the ellipticity e and the cosine of the angle 
between k and B q . The calculation imposes the requirement that k has a 
positive projection in the +R direction. The normalized magnetic helicity and 
the degree of polarization track each other very well. Also note that the 
ellipticity is large when both the magnetic helicity and the degree of 
polarization are large, implying nearly circular polarization. Therefore, 
subject to the condition that the fluctuations have phase speed less than the 
solar wind speed, the waves must be left-hand circularly polarized in the 
plasma frame [Smith et al. 1983]. In addition, the cosine of the angle 
between k and B q , cose , is large for frequencies from 0.001 to 0.1 Hz. Thus 
these fluctuations should not be compressive. Because the rest frame polariza- 
tion is left-hand, these fluctuations are probably AlfvSn waves propagating 
nearly parallel (or anti-parallel) to the ambient magnetic field. 

We have checked that these waves are noncompressive by correlating the 
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density measurements p with the magnetic field magnitude |b|. Alfven waves 
(and parallel propagating fast mode waves) should show little if any correla- 
tlon between | B l and o if 4 * 0. This correlation is shown in Figure 4b using 
14 degrees of freedom for frequencies from 0.025 to 0.04 Hz, which is the only 
range that overlaps the magnetic spectrum. Note that the magnetic field 
magnitude and plasma density correlation appears to give almost no correlation 
since the peaks in this frequency range tend to oscillate about zero. However 
this should be interpreted with some caution because the fluctuations have 
very small amplitude and the correlation may not be well resolved. 

The eigenvalue analysis can only determine the wave phase velocity 
direction to within a sign. However, this ambiguity can be resolved by 
calculating the cros s-helicity spectrum. The reduced normalized cross- 
helicity spectrum is defined as <? c (k) “ 2H c (k)/E(k) where E(k) and K c (k) are 
the spectral decomposition of the total energy spectrum (magnetic plus 
kinetic) and the cross-helicity spectrum, respectively. The cross helicity 
measures the correlation between the velocity and magnetic fluctuations and is 
defined by: 


H » /d s x v*b * /d 3 k H (k) 
c c 

where b “ <5B//(4irp). The total energy in these Alfven velocity units is 

E- J“d 3 x (v 2 + b 2 ) * /d 3 k [E v (k) + E b (k)] 

When a c (k) is near ±1, Alfvenic fluctuations are present in the data. Thus, 
if the magnetic field and the velocity fluctuations are in the same direction 
(a = +1) the wave energy is propagating antiparallel to the mean magnetic 
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field. If they are in the opposite sense (a - -1), then the wave energy is 

c 

propagating parallel to Figure 4a shows the normalized cross-helicity 
o c (k) spectrum as a function of frequency. The spfifttrum is negative and quite 
large (o c « -0.95) for frequencies in the range of 0.025 to 0.04 Hz. The 
negative sign indicates that the fluctuations are propagating parallel to the 
average magnetic field which is directed away of the shock and toward the Sun. 
The fact that the cross-helicity is large also tends to confirm that the 
fluctuations are Alfvenic and not whistler waves because the plasma data used 
only included protons. 

Further confirmation that these fluctuations are MHD and not whistler 
waves is indicated by the Alfven ratio r A (k) - E v (k)/E fe (k) plotted in Figure 
4c. This quantity measures the degree of equipartition between the magnetic 
energy and the kinetic energy. Kraichnan £1965] predicted that in fully 
developed MHD turbulence the two should be approximately equal at wave numbers 
large compared to those characterizing the energy-containing scales. Note 
that for whistler waves, r^ should be nearly zero because the plasma measure- 
ments we use do not include the electron data [Matthaeus and Goldstein, 1982; 
Goldstein et al., 1983]. 

January 29 — Downstream 

A similar analysis has been performed in the downstream interval (from 
0800:45 to 0914:39 UT) of the January 29 quasi-parallel reverse shock. In 
this interval the components of the average magnetic field are (-1.47, -0.76, 
-1.07) nT and the average solar wind speed is 337 km/s. The typical amplitude 
of these fluctuations is again sizeable, ( <$B) mg / 1 <B q > I = 0.42. The magnetic 
field power spectrum (Figure 5) has an f * power-law dependence at high freq- 
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uencles and It contains a peak centered about f 3 0-013 Hz indicating the 
presence of quasi-monochromatic fluctuations in the downstream region. 

Figures 6a-d show the calculations of the magnetic heliclty spectrum <x^, 

the degree of polarization D, the ellipticity e and the cosine '>f the angle 9 

between k and B - Note that near f 3 0.013 Hz the magnetic heliclty is 
o 

positive (a =• 0.3) indicating left-hand helices. Similarly the degree of 
m 

polarization and ellipticity yield D 3 0.6 and e « 0.5 respectively. The 
fluctuations near f 3 0.013 Hz appear to be elliptically polarized Alfven 
waves propagating parallel to the magnetic field as indicated In Figure 6d. 

In Figure 7a-c the calculation of the cross-helicity spectrum, the |B|-p 
correlation and the Alfven ratio r^ are shown « Because the cross-helicity is 
negative and large (a c 3 -0.8), the fluctuations are propagating parallel to 
the magnetic field (which in this case is toward the shock and the Sun) . The 
magnitude of the cross-helicity again suggests that these are Alfvenic 
fluctuations and not whistlers. Because the rest frame polarization is left- 
hand (c > 0) , we infer that the fluctuations are Alfven modes propagating 
quasi-parallel to the ambient magnetic field. Figure 7b shows the | B I — P 
correlation using 14 degrees of freedom for frequencies near f 3 0.013 Hz 
which are the only range that overlaps with the magnetic spectrum. Note that 
this calculation appears to indicate almost no correlation since the peaks in 
this frequency range tend to fluctuate about zero. The Alfven ratio r^ for 
this interval is plotted in Figure 7c. At all frequencies r A is near one as 
expected for MHD fluctuations. 

Although the peak at f - 0.013 Hz is the most prominant (see Figure 5), 
there is also a relatively broad enhancement at higher frequencies . From 
Figures 6a-d we see that becomes negative (right hand helices) at higher 

frequencies. The maximimum value is reached at about f 3 0.14 Hz where a 3 
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-0,6. The good agreement betwiu&n a m , D and s for frequencies about 0,05 - 0.2 
Hz suggests that the downstream region may contain both Alfven and fast mode 
waves. However we are unable to determine with certainty if these right-hand 
waves are indeed fast MHD modes because the Nyquist frequency of the plasma 
data is well below the frequency range of these waves . 

February 3, 1978 

Another example of low frequency waves upstream of an interplanetary quasi- 
parallel forward shock is presented in Figure 8. Here we show a forward shock 
forming at about 1928 UT on February 3, 1978, We have selected the upstream 
region (excluding the shock) which corresponds to the subinterval from 1850 to 
1926 UT for the spectral analysis. During this subinterval the components of 
the mean magnetic field are (-1.60, 1.56, -1.16) nT and the average solar wind 
speed is 399 km/s. The fluctuations in this example have about the same 
magnitude as before; ( <SB ) ritts /l < B 0 >I “ 0, 26. The shock normal components for 
this event are n » (-0.69, 0.62, -0.366) and the angle between the normal and 
the average upstream magnetic field is about 7.57°. 

In Figure 9 we show the magnetic power spectrum for this event which has 
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an f ' power- law dependence at high frequencies . A bulge in the power 
spectrum centered about f » 0.035 Hz indicates the presence of quasi-monochro- 
matic fluctuations. a_ (k), D, e, and cos9 are plotted in Figure lOa-d. The 
magnetic hellcity decreases with frequency becoming negative (g^ a -0.5) in 
the range 0.025 - 0.17 Hz, indicating the presence of right-hand helices. 
Also, D and cr m (k) track each other well. The ellipticity as shown in Figure 
10c is also large (e = 0.9). Thus, these fluctuations are nearly circularly 
polarized . 


-12 


Because the rest frame polarization Indicates right-hand polarization, 
these fluctuations are probably fast MHD modes propagating almost exactly 
parallel to the ambient magnetic field (cf. Figure lOd) . To determine the 
sign of the direction of propagation, we evaluated the cross helicity (Figure 
ila) . From the positive values in this frequency range, we conclude that 
these fluctuations are propagating away from the shock (and the Sun). Protons 
streaming away from the sujUc can thus be in resonance with these waves. The 
large values (r » 1) of the Alfven ratio in Figure 11c lend further con- 
firmation that these fluctuations are MHD and not whistler. 

However, there is an Interesting difficulty with this interpretation in 
terms of linear wave modes. Note that the |b|~p correlation (Figure lib) is 
relatively large (=* 0-6). This is rather surprising because if 9 is very 
small, the MJS'-Sn and the fast MHD branches become degenerate to lowest order, 
and <( on e‘t& basis of linear theory, little compression should be present. 

We have investigated this phenomenon more fully. In Figures 12a and 12b, 
the power spectrum of the density fluctuations and the power in { B | are 
plotted. Both spectra show an enhancement between 0.025 and 0.04 Hz, confirm- 
ing the high correlation noted in Fig. 11b. However, the amount of power in 
the spectrum of f B ( is much lower than in the components (Fig. 9). The 
observed correlation between I B I and P, therefore, may be a rather high order 
effect . 

It is well known that the degeneracy of the Alfven mode and the fast mode 
is broken if one treats the fast mode to second order in wave amplitude. 
Barnes and Hollweg [1974] studied an aspect of this problem. If one assumes 
that to lowest order the MHD waves have linear polarization, then to second 
order the relationship between density fluctuations and transverse magnetic 
field fluctuations for nearly parallel propagation is given by: 
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«0 («B a ) 2 - <6B A 2 > 

_ . _i : — (1) 

o 2(B * - 4iry P ) 

0 o p 0 

where <5p is the density perturbation in the wave, p Q is the mean density, <5B A 
is the amplitude of the transverse magnetic wave field, <> represents a 
temporal average, and P Q is the mean pressure. 

The difficulty with applying (1) to our observations is that the observed 
wave is nearly purely circularly polarized. For circularly polarized fluctua- 
tions, there is no second order correction for parallel propagation and the 
right hand side of (1) vanishes. Another related interpretation of the*, data 
is that the jBj-p correlation represents not a high order correction to the 
nHB mode amplitude,, but rather evidence for mode coupling between the MHD wave 
thought of as a pump and another ‘'daughter” wave via a raodulational instabi- 
lity. The theory of modulational instabilities has been developed by Lashmore- 
Dsvias [1976], Derby [1978], Goldstein [1978], Sakai and Sonnerup [1983], 
among many others. If one treats the pump wave as being circularly polarized, 
and assumes parallel propagation, then the wave equation relating density to 
magnetic field is 

3 2 <5o _ 2 3 1 6p * 1 3 2 (H t »<$Bj) , _JL_ 3 2 (6B,»6B,) 

at 2 s 3 z 2 4ir 3 z 2 8ir gz 2 s "'’ 

If we now assume that the daughter wave is also linearly polarized, the 

last term on the right-hand-side of (2) is zero. An analysis similar to that 
which led to (1) shows that the relationship between <5p, Bj_ and <sBj_ is 
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6p (2 B x 63 x ») 

a . £ ^ ^ 

p o M 0 * - ^Y p P 0 ) 

We have estimated the extent to which (3) is satisfied for this data 
interval. The, quantity 5p is estimated from th/^ power spectrum of the density 
by integrating the power in the frequency range 0.025 - 0.04 Hz. One can 
estimate (Bj_ SB^) by rotating the magnetic field data into the mean-field 
coordinate system. The components transverse to the mean are then easily 
extracted and SBj_(t) * [B x (t) - <Bj>] . 

We take y * 5/3, P * NKT and T * 8.26 x 10* °K, which allows both slde3 
p ’ o ’ 

_2 

of (3) to be estimated independently. The result is that <Sp /p q « 4.7 x 10 , 

-3 

while the right-hand- side of (3) is 8.4 x 10 . The fact that the two differ 

by a factor of six suggests that either additional physical processes are 
important or that our estimates of the quantities in (3) are inadequante. The 
latter is certainly possible in that, first of all, the interacting waves are 
not monochromatic as was assumed in (3); although the power spectrum is peaked 
' near 0.035Hz, it is certainly not a delta function (see Figure 12). Secondly, 
the assumption that the magnetic and density fluctuations are linearly 
polarized can be only approximately valid. Finally, (3) was derived under the 
assumption that the pressure is isotropic, but the presence of pressure 
anisotropies may significantly modify the analysis (see Hollweg, 1971). 

3. Linear Instability Analysis 

At planetary shocks similar waves as those described in the previous 
section have been investigated by Barnes (1970'], Gary at al. (19811, Sentman 
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et al. [1981J, Lee (1982], and Goldstein et al. [1983], among others. The 
principal result of these papers is that such waves may be generated by either 
resonant particle interactions excited by the protons reflected (or accelerat- 
ed) at' planetary shocks, or by a iion-resonant interaction of the waves with a 
"diffuse" ion population which airlses from pitch-angle scattering of the 
reflected component. In this section, we present a similar analysis and 
discussion of the generation mechanism for our observations. 

In this paper we limit our discussion to the field aligned electromagnetic 
ion beam instability. A more complete analysis, including oblique propagation 
and a comprehensive search of the Voyager plasma data for evidence of the 
existence of the required ion distributions will be deferred to a later paper. 
Our model consists of a background electron-proton plasma that is uniform and 
infinite through which an ion beam is streaming. We assume that the frequency 
o of the fluctuations (in the plasma frame) is smaller than the ion gyrofreq- 
uency The drift velocities of both the background ion and the ion beam 
(V q £ and V ob , respectively) are also directed along the magnetic field, but 
the electron drift velocity is zero (V Qe ■ 0). The medium is assumed to be 
charge neutral (n fi - n^+ n^) and to have zero current (n i V Qb * 0) 
where n^, n^ and n^ are the number densities of the background electrons, ions 
and the beam components, respectively. The zeroth order distribution for the 
jth component is a drifting bi-Maxwellian of the form 



(4) 


where and a„j are the perpendicular and parallel thermal velocities 

defined as 


J 
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kg Is Che Boltzmann constant, is t'he mass of the jth species and T,,^ and 
T x ^ are the parallel and perpendicular temperatures of the jth species 
respectively. 

The dispersion relation for parallel propagating low frequency electro- 
magnetic waves is [see, for example, Montgomery and Tidman (1964)]: 


1 2 2 2 2 Vj - X 1 (v, ” v x> 

* •- u - k„ c + irZ u / S dv„ f i dv, - 0 

3 pj 


<j-k„v„±n 


j 


(5) 


and 


Xj( v "» v j,)" (u-k„v„> 


3f 


(o) 


3f 


( 0 ) 


+ k..v, 


3v. 


3v„ 


2 1/2 

where u is the (complex) wave frequency, * (4™^ /up is the plasma 
frequency of the jth species, c is the velocity of light and 52 ^ * q^Bpinp. 
The plus and minus sign in the denominator of (5) refers to right and left 
hand polarization, respectively. After substitution of (4) into (5) we obtain 
the dispersion relation for the right and left hand polarized waves [see, for 
example, Scharer and Trievelpiece , 1967; Sentman et al., 1.981]: 


D. 


R u 2 - k„ 2 c 2 + S u pj 2 Upp) + Ap Up/21 - 0 


u-k„V 




oj 


52 . 


= U ± 


k„a„ 


j 


k„ ct, t ^ 


( 6 ) 


Aj = l-c i //a„j is the anisotropy, and Z(S) and Z'U) are the plasma disper- 
sion function and its derivative with respect to 5 [Fried and Conte, 1961]. 
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To apply this formalism to the observations, we first estimate the plasma 
frame frequency and wave number using the Doppler shift equation together with 
the assumption that u « The eigenvalue analysis provides both 9 and the 
angle between £ and $ » Because the rest frame polarization is known from 
the magnetic helicity spectrum, there are only two possibilities to consider, 
viz., whether the instability is resonant or nonresonant. We have investigat- 
ed both possibilities for each interval. 

The overall results of this instability analysis can be summarized in a 
general way before presenting th& detailed computations below. The cold 
plasma, resonant Instability (u - kwV^ + ^ " 0) was recently discussed by 
Goldstein et al. [1983]. For parallel propagation, there is only one unstable 
mode, viz., the right-hand fast mode, in which the beam velocity and wave 
phase velocity are parallel. In addition, there is a nonresonant instability, 
also right-hand, for which the beam and wave propagate antiparallel [see, 
e.g., Sentman, et al., 1981; Gary et al., 1984]. When the background plasma 
has a finite temperature, two new instabilities appear that are left hand 
polarized. First, there is a resonant instability (<j - h„V o j j - 0 « 0) that 
grows only if A, > 1. In this case the beam velocity and wave phase velocity 

D 

are antiparallel. Finally, there is a nonresonant left-hand instability (with 
the beam and wave propagating parallel) that can grow if the beam temperature 
is sufficiently large. can be zero for this mode. 

January 29, 1978 — Upstream 

The fluctuations in the upstream region of the January 29 shock were left- 
hand polarized in the range 0.025 - 0.07 Hz and were propagating parallel to 
the magnetic field, away from the reverse shock (into the upstream region). 


f 
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The waves must be Doppler shifted "down" to u 1 * u - &*V where o r is the 
wave frequency as measured in the spacecraft frame in radians/s. By combining 
this with (6), we can determine the frequency of the wave in the solar wind 
frame, its phase velocity, its wavelength and the resonant velocity of the 
particles for our observations. 

The plasma parameters for the background plasma, estimated from the 
Voyager data, are summarized in Table 1. Unstable roots of (6) were found 
under the assumption that the Interaction was either resonant or non-resonant . 
The beam parameters used in the calculation of the wave characteristics were 
determined by requiring both that u < ftp and that the wavelengths and freq- 
uencies were in the observed range with the maximum growth rate occuring 
inside this range. 

First consider the resonant interaction. Note that for this to be the 
physical mechanism for excitation of the observed waves, the beam will have to 
be propagating toward the shock front. Thus, we expect that this is not a 
likely scenario, but it is instructive to explore the range of physical 

parameters required to excite this mode. The unstable waves are found to have 

—3 

frequencies (2.9 - 6.3) * 10 Hz and wavelengths between (4.82 - 13.7) * 10 s 

cm. The maximum normalized growth rate (y/D ) for the resonant interaction is 

P 

—2 —2 

2.34 x 10 at a wavelength of 7.13 * 10* cm and u ■ 3.09 * 10 rad/s. The 

-2 

average gyrofrequency for this event is 2.04 x 10 Hz (Table 1). Note, from 
Table 1, that this instability requires a very large ion beam thermal anisotro- 
py (=5) and a beam velocity of V ^ » -41.25 km/s (toward the shock). 

In the case of the nonresonant interaction the unstable waves and beam 
will both be propagating away from the shock, which would appear to be a more 

reasonable situation than the geometry of the resonant interaction. The 

-3 

unstable waves are in the range of (2.7 - 6.2) x 10 Hz with similar wave- 


lengths as for the resonant interaction. The maximum normalized growth rate 

—3 8 

for the non-resonant interaction is 2.68 x 10 at a wavelength of 6.56 * 10 

—2 

cm and frequency u » 3.2 x 10 rad/s. For this instability, a hot isotropic 
"diffuse” proton population of about 222 ev streaming away of the shock at a 
velocity of V ^ * 103 km/s is required, which is just barely sufficient for 
the beam to propagate back upstream. This instability is insensitive to the 
value of A^. 

Recently, Hoppe and Russell [1983] have reported evidence for left-hand 
Alfven wave fluctuations in the earth's foreshock region from ISEE-1 and 2 
data. They concluded that the Alfven mode could have been excited only by a 
diffuse population and not by a resonant beam interaction. Goldstein at al. 
[1983] argued that this was not necessary, because for oblique propagation the 
left-hand Alfven mode can be resonantly excited. The situation that we have 
here is complementary in that we find that with a sufficient beam anisotropy 
amplification of left hand Alfven waves is possible as anticipated by Hoppe 
and Russell [1983]. 


January 29, 1978 — Downstream 

A similar analysis has been performed in the downstream interval of the 
January 29 event for the observed frequency range of 0.01 - 0.015 Hz. During 
this interval left-hand waves are propagating parallel to the magnetic field 
toward the reverse shock. Thus, the frequency in the spacecraft frame must be 
Doppler shifted "down" by u' * u - &*V gw . Table 2 shows some of the results 
obtained assuming either a resonant and a nonresonant interaction. 

For the resonant interaction the unstable modes are in the range of (1.1 - 

-3 q 

2.0) x 10 Hz with wavelengths in the range (1.35 - 2.7) x 10 cm. The 

-4 9 

maximum normalized growth rate is 5.71 x 10 at a wavelength of 1.8 x 10 cm 
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_3 

and frequency o * 9.45 x 10 rad/s. The beam parameters used in this 
calculation were ■ 1.7 and V ^ * -184 km/s. 

Similarly, the unstable modes for the non-resonant interaction occurs in 

_3 

the range of <1.3 - 2.4) x 10 Hz for the same wavelengths as in the resonant 

-3 

case. The maximum normalized growth rate for this interaction is 5.36 x 10 

—3 

at a wavelength of 1.67 x 10* cm and frequency u - 2.0 x 10 rad/s. For this 
particular event the resonant interaction is much more probable because the 
ions are streaming away from the shock (into the downstream region) with 
modest anisotropies that are actually typical of the ambient solar wind. 

At higher frequencies (Figure 6a) , (0.05 - 0.2) Hz we found field aligned 
right-hand polarized waves. In the absence of cross helicity data in this 
frequency range, we shall assume that we are dealing with fast mode MHD waves 
propagating toward the shock (and the Sun). The dispersion relation is then 
given by (6) using the plus sign for right-hand waves. Since the solar wind 
is convecting away from the Sun, the waves observed at the spacecraft position 
must again be Doppler shifted "down'' . The instability analysis indicates that 
in the plasma frame u extended from frequencies below the ion gyrofrequency to 
frequencies above it. This suggests that these waves may be right-hand ion 
cyclotron waves propagating nearly parallel to the magnetic field. These 
modes have been observed near interplanetary shocks by Tsurutani et al. 
[1983] . 

February 3, 1978 

In this case the fluctuations were right-hand with frequencies in the 
range (0.025 - 0.17) Hz (in the spacecraft frame) propagating nearly anti- 
parallel to the magnetic field away from the forward shock (into the upstream 
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region) . Similarly the wave characteristics in the plasma frame have also 

been calculated, using the dispersion relation (6) (with the plus sign for 

right hand waves) and the Doppler shift w’ « u + k* v sw * Table 3 contains some 

results for the beam parameters that best fit the observations for both the 

resonant and nonresonant wave-particle interactions. In the resonant case, 

-3 

the wave frequencies in the plasma frame are about (4.4 - 26) * 10 Hz, with 

wavelengths between (2.19 - 9.51) x 10* cm. The maximum normalized growth 

-2 

rate (y/Q^) for resonant interaction ocurred at 3.49 x 10 at a wavelength of 

—2 

3.42 x 10 8 cm and frequency u - 9.35 x 10 rad/s. These modes can be excited 
by a proton beam streaming away from the shock (into the upstream region) with 
a thermal spread of about 38 ev and beam speed of -150 km/ s (Table 3). 

Excitation of the nonresonant mode requires an energetic proton beam 
propagating back toward the shock (and the sun) . This is a very unlikely 
situation, and we futher find that a very high beam speed of 557 km/ s is 
required. The frequencies in the plasma frame of these unstable modes are 
between (1.65 - 13.0) x 10 Hz for the same wavelengths obtained in the 

resonant case and the maximum normalized growth rate for this case is 2.62 x 

•••l. -2 

10 at a wavelength of 3.57 x 10* cm and u * 4.0 x 10 rad/s« 

4. Discussion and Summary 

We have presented three examples of low frequency waves associated with 
interplanetary shocks; a quasi-parallel forward shock and a reverse shocks. 
Two MHD modes have been identified. The first is a fast mode with character- 
istic frequencies in the range of 0.025 - 0.17 Hz (in the spacecraft frame). 
This mode is right-hand elliptically polarized propagating along the magnetic 
field. The second is an Alfv&n mode with characteristic frequencies in the 
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range of 0.025 - 0.07 Hz (in the spacecraft frame). This mode is left-hand 
polarized and propagates within a few degrees of the magnetic field. The 
magnetic power spectra observed in the vicinity of these shocks are much 
steeper at high frequencies than is characteristic of either the ambient solar 
wind and or planetary shocks . 

The analysis we have presented describes some of the basic properties of 
waves in the vicinity of interplanetary shocks. Our observations are consis- 
tent with either a resonant or nonresonant electromagnetic ion beam instabili- 
ty. In two of the three situations, the February 3, 1978, and the downstream 
side of the January 29, 1978 shocks, the geometry may be favorable for 
detectirs the reflected proton distributions in the plasma data. This effort 
will be reported in a subsequent paper. 

In spite of the limitations of the linear theory, and in spite of the many 
assumptions and limitations of our theoretical treatment, the electromagnetic 
ion beam instability appears capable of accounting for the major features of 
the observations. Similar conclusions have been reached in connection with 
planetary flow shocks by Barnes [1970], Gary et al. [1982], Lee [1982, 1983] 
and Goldstein et al. [1983], among others, and by Tsurutani et al. [1983] in 
their analysis of interplanetary shocks . 
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Figure Captions 

Figure 1. Voyager 2 magnwtic field and plasma data of a reverse quasi- 
parallel shock occuring at 0918:17 UT on January 29, 1978 at about 2.16 AU. 

Figure 2. Trace of the magnetic energy spectrum from the upstream region 
(0915k 10-1000 UT) fot the January 29 event. The data used are 1.92 s averages 
and the spectrum contains 26 degrees of freedom. 

Figure 3. a. The normalized reduced magnetic helicity spectrum after 
Matthaeus and Goldstein [1982J. 

b. The degree of polarization D as obtained from an eigenvalue 
analysis of this interval. 

c. The ellipticity e as calculated from the eigenvalue analysis. 

d. The cosine of the angle 9 between k and B Q as obtained from 
the eigenvalue analysis imposing the requi ment that k has a positive 
projection in the +R direction. 

Figure 4. a. The normalized cross-helieity spectrum a after Matthaeus and 
Goldstein [1982]. The data used are 12 s averages and the spectrum contains 
14 degrees of freedom. 

b. The correlation between I B| and p. 

c. The Alfven ration r A - E v (k)/E^(k) where E^Ck) and E fe (k) are 
the reduced spectra of kinetic and magnetic energy respectively. 
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Figure 5 . The trace of the magnetic power spectral density for the 
downstream interval (0800:45-0914:39 UT) on January 29 event. This spectrum 
has 26 degrees of freedom. 

Figure 6. Plotted in a-d are o , D, e, and cosO, respectively, for the 
downstream region. 

Figure 7. Plotted in a-c are a the | B ! — p correlation, and the Alfven 
ratio r^ respectively, for the downstream region. 

Figure 8 . Voyager 1 magnetic field and plasma data of a quasi-parallel 
forward shock (in formation) occuring at 1928 UT on February 3, 1978 at about 
2.25 AU. 

Figure 9. Trace of the magnetic energy spectrum for the upstream region 
(1850-1926 UT) of the February 3 event. 

Figure 10. The same as Fig. 6, but for the upstream region of the February 
3 event. 

Figure 11. The same as Fig. 7, but for the upstream region of the February 
3 event . 

Figure 12. a. The density spectrum for the February 3 event containing 14 
degrees of freeedom. 

b. The magnetic field magnitude spectrum. 


Table 1. Plasma parameters for January 29, 1978 (Upstream) 
Left-hand, Nonresonant Instability 


Background electrons 

Background ions 

Ion 

beam 

Resonant 

Nonresonant 

Resonant 

Nonresonant 

Resonant 

Nonresonant 

n^ (///cm*) .50 

.50 

0.47 

0.47 

O'. 005 

0.005 

k B T Mj (eV) 6.15 

6.15 

1.47 

1.47 

26.65 

222. 

kjjTj. j (eV) 6.15 

6.15 

1.47 

1.47 

133.26 

222. 

V a (tan/ a) 0.0 

oj J 

0.0 

0.42 

-1.04 

-42.25 

103. 

The Alfven speed is 4. 

12 x 10*, V 

’ sw 

* 3.66 x 

10 7 , 0*2. 
P 

04 x 10~ 2 



Table 

2. Plasma parameters for January 29, 1978 
Left-hand, Resonant Instability 

Background electrons Background ions 

(Downstream) 
Ion beam 


Resonant Nonresonant 

Resonant Nonresonant 

Resonant 

Nonresonant 

n^ (///cm 3 ) 

1.36 

1.36 

1.27 

1.27 

0.014 

0.014 

k B T "j (eV) 

6.09 

6.09 

5.52 

5.52 

141.84 

319.14 

(feV) 

6.09 

6.09 

5.52 

5.52 

248.22 

319.14 

(tan/ a) 

0.0 

0,0 

1.86 

0.93 

-184.3 

92.15 

The Alfven speed is 3.68 x 

10 s , V 
* sw 

1 3.37 x 

10 7 , a = 3.1 
* P 

01 x 10~ 2 



. // 
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Table 3. Plasma parameters for Feburary 3, 1978 (Upstream) 
Right-hand , Resonant Instability 


Background electrons Background ions Ion beam 



Resonant 

Nonresonant 

Resonant 

Nonresonant 

Resonant 

Nonre/sonant 

(///cm’) 

1.64 

1.64 

1.52 

1.52 

0.016 

0.082 

k^T„j (eV) 

5.27 

5.27 

5.64 

5.64 

38.4 

9.6 

V*} (eV) 

5.27 

5.27 

5.64 

5.64 

38.4 

9.6 

V oj (km/s) 

0.0 

0.0 

1.52 

-29.34 

-150, 

557.4 

The Alfven speed is 4.28 x 10*, V 

sw 

- 3.99 x 

10 ? , Q * 3. 
P 

84 x 10“ 2 
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